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circular polarisation components, have been impeded due to non-optimised 
instrumentation. The need for suitable observing conditions and the availabil¬ 
ity of luminous targets are also limiting factors. The science motivation of any 
instrument adds constraints to its operation such as high signal-to-noise (SNR) 
and detector readout speeds. These factors in particular lead to a wide range 
of sources that have yet to be observed. The Galway Astronomical Stokes Po- 
larimeter (GASP) has been specifically designed to make observations of these 
sources. 

GASP uses division of amplitude polarimeter (DOAP) ([Compain and Dr^ 


villon 19981 to measure the four components of the Stokes vector (I, Q, U and 


V) simultaneously, which eliminates the constraints placed upon the need for 
moving parts during observation, and offers a real-time complete measurement 
of polarisation. Results from the GASP calibration are presented in this work 
for both a ID detector system, and a pixel-by-pixel analysis on a 2D detector 
system. 


Following Gompain et al. (1999) we use the Eigenvalue Galibration Method 


(EGM) to measure the polarimetric limitations of the instrument for each of 
the two systems. Gonsequently, the EGM is able to compensate for systematic 
errors introduced by the calibration optics, and it also accounts for all opti¬ 
cal elements of the polarimeter in the output. Initial laboratory results of the 
EGM are presented, using APD detectors, where errors of 0.2% and 0.1° were 
measured for the degree of linear polarisation (DOLP) and polarisation angle 
(PA) respectively. Channel-to-channel image registration is an important as¬ 
pect of 2-D polarimetry. We present our calibration results of the measured 
Mueller matrix of each sample, used by the EGM, when 2 Andor iXon Ultra 
897 detectors were loaned to the project. A set of Zenith flat-held images were 
recorded during an observing campaign at the Palomar 200 inch telescope in 
November 2012. From these we show the polarimetric errors from the spatial 
polarimetry indicating both the stability and absolute accuracy of GASP. 


Keywords polarisation • EGM • calibration • APD • EMGGD. 


1 Introduction 

A Stokes polarimeter is described as ‘complete’ when the polarisation state 
analyser (PSA) can measure the four linearly independent states of polari¬ 
sation of light. The calibration of a polarimetric tool is a two-step process. 
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Firstly, the optical elements of the system must be precisely orientated with 
very careful alignment, and secondly, the behaviour of these optical elements 


must be easily described and characterised by physical properties (Compain 


et al. 19991. Therefore, the combination of several measurements in different 


configurations, to determine the main values and shortcomings. Any optical 
element in the light path that has not been fully characterised, can introduce 
modifications to the state of polarisation of light and bias the estimation of 
the Stokes parameters. 

The linear component of optical polarised light emission has been well 
documented, however, observations have been limited due to instrumental 
considerations to periods of excellent observing conditions, and to steady (or 
slowly or periodically-varying) sources (Collins et al., 20131. In astronomical 
polarimetry, the optical circular component is less established. One important 
consideration in the development of any instrument, particularly in this work, 
is the ability to achieve adequate SNR for the science action, or minimise all 
sources of instrumental error. 

GASP was designed as a DOAP polarimeter using a modified Fresnel 


rhomb (Compain and Drevillon, 19981. Collins et al. (20131 reported an initial 
design of GASP and results based on a pseudo-inverse calibration. In order to 
calibrate a polarimeter, the polarimetric optics/materials calibrating the po- 
larimeters must also be calibrated. The orientation of the elements that make 
up the polarisation state generator (PSG) - the angle of transmission axis of 
the polariser and the fast axis of the waveplate - must be calculated. If there 
are any errors in the PSG, used to calibrate/train the PSA, these errors will 
also propagate through to the instrument system matrix. 

The ECM has the ability to calibrate a Stokes polarimeter in the exact 
configuration in which it will be used. It is also effective in the case where the 
PSG is a required output of the calibration process, which can be found in 


a number of references including Compain et al. (19991; de Martino (20041. 
Another advantage of this method is that it uses simple, off the shelf glass 
optics that are characterised, and calibrated, during the ECM. The ECM can 
independently calibrate the PSA and PSG, and works in such a way that the 
instrument matrix can never be polluted by errors (described above) from 
the PSG. The main science motivation for this instrument is time resolved 
observations of the Crab pulsar, which has a period of approximately 33 ms and 
magnitude of only 16.5. These characteristics require the removal of all/any 
instrumental error, including that of calibration. 
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We report the operation of GASP using an ECM approach, using two differ¬ 
ent optical-detector systems. The behaviour of any optical or electronic device 
can be distorted by first-order failures, as in the case of a photo-elastic mod¬ 
ulator (Drevillon 19931. Therefore, several measurements must be combined 
in a number of different configurations to determine the main values and the 
shortcomings. The optical set-up is then adapted to account for these measure¬ 
ments. These aspects of the calibration process have already been reviewed in 
the literature (Hauge 1978 Thompson et al., 1980 Azzam and Lopez 19891. 
Furthermore, one can notice that usual calibration procedures are difficult to 
implement in-situ because of the influence of the optical elements that are 
necessarily included in the light path (filters, windows, lenses, mirrors, etc.) 
All these elements can induce modifications of the state of polarisation of light. 


It has been discussed that the ECM minimises the main sources of er¬ 
ror generally associated with calibration where it requires a complete system 
(the polarimeter), a Polarisation State Generator (PSG) and 4 samples. No 
assumptions are made to the system, the PSG and Mueller matrices of each 
sample are actually characterised by the ECM. The system matrix is simply a 
function of the wavelength and must be consistent from calibration to science. 


The general form for the matrices A and W are n x 4 and 4 x m, respec¬ 
tively, where n and m are greater than 4. n > 4 mean that more than four 
input Stokes vectors are generated, and m > 4 means that the output Stokes 
vectors are over-determined. This is the case for GASP. The ECM is fully 
compatible with other size matrices and the dimensions of the intermediate 


matrices must be adapted (Compain et al. 19991. 


An investigation of the ECM, using the two configurations of GASP, one 
using Avalanche Photodiodes (APDs), and one using imaging detectors, will 
be presented in two stages. First, the APDs were used to record data from 
a laboratory experiment. Then 2D ECM data was obtained during a 4 night 
observing run at the 200 inch Hale telescope at Palomar Observatory. The 
project was loaned 2 Andor iXon Ultra 897 EMCCD (Electron Multiplying 
Charge Coupled Device) detectors for the 4 night campaign. As well as an 
experimental calibration run, GASP was also used to observe the Crab pulsar 
(among other targets) at a frame rate > 1000 Hz to obtain a phase-resolved 
polarimetric signal. 
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2 The Galway Astronomical Stokes Polarimeter 


GASP uses a modified DOAP design, which was implemented in the laboratory 
by Compain and Drevillon (19981. A full description of the polarimeter, in 
its initial design stage, is found in Collins et al.| ( |2013 ), with further optical 
development detailed in Kyne ( 2014[ ). Figuregives a general description of 
the polarimeter 
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Fig. 1 The optical layout for a 2D detector system using GASP; this is a general design 
that can be adapted for a specific telescope. The reflected path (RP) is found on the top 
right hand corner, and the transmitted path (TP) is found at the bottom right hand corner. 
The optical prescription (with the exception of the main beam splitting prism) changes 
depending on the choice of telescope and detector. 


2.1 The Polarisation State Analyser 

The Polarisation State Analyser, PSA, is the instrument used to estimate the 
Stokes vector of a beam of light. In the case of an imaging polarimeter, the 
PSA is used to estimate a separate Stokes vector for each field position in the 
image. In GASP, the PSA is the combination of the main beamsplitting Fresnel 
rhomb, two polarising beamsplitters such as Wollaston or Foster prisms, the 
four simultaneous detectors, and any optical elements between the prism and 
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the detectors. The PSA in GASP was originally designed as an 8-detector/lD 
system, using Foster prisms, where FOV size and image quality were not of 
great importance. The Wollaston (Foster) prisms must be aligned to an integer 
multiple of 45°, with respect to the optical axis of the instrument. A 45° angle 
produces an invertible calibration matrix and splits each beam 50:50, once the 
system has been aligned experimentally; it must be calibrated based on this 
alignment. The PSA, with its base in the plane of incidence, acts as a first- 
stage beamsplitter and provides the essential polarimetric properties. These 
properties are displayed, as usual, by the reflection and transmission Fresnel 
coefficients of the prism (Azzam and Bashara 19881, 


The PSA has a characteristic matrix. A, that fully describes the instrument 
polarimetrically, i.e., it contains all the information for how well the instrument 
is calibrated, experimentally. This matrix is the output of the ECM and con¬ 
tains both the system and demodulation matrix, which can be described best 
as a theoretically corrected system matrix, which is a better representation of 
the hardware used in constructing GASP. Another inclusion to the A matrix is 
the detector gain; the EGM allows for a different gain on each detector/pixel. 
The final Stokes vector S, is given by the following relationship: 


S = A-^I 


( 1 ) 


where A is the 4x4 experimental system matrix and I is the 4 x n intensity 
array, n can be of any length depending on the acquisition time of the data. S 
will be the same size as I. 


Collins et al. (20091 describes a design where the four components of the 


Stokes vector are measured simultaneously with a retardance error ^1%, over 
a spectral range of 400 - 800 nm. The DOAP configuration contains no moving 
or modulating components, where it is possible to achieve high time resolution 
(of order microseconds), and temporal stability based on the choice of detec¬ 
tor. The optical configuration of GASP can vary for reasons explained in |Kyne| 
( 2014[ ), depending on instrument size and detector availability. The values in 
Table are the properties corresponding to the current GASP rhomb and con¬ 
figuration for the instrument which observed at Palomar. All measurements, 
theoretical and experimental, use these values. 
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Table 1 Geometrical and polarimetric properties of the current GASP prism-rhomb from 
|Collins et aL|||2013|l and|Kyne|l|2014||. 

Prism Geometry 

Glass Properties 

Polarimetric Properties 



Reflection 

Transmission 

01 = 78.50° 

A = 589nm 

R = 0.3533 

T = 0.4352 

02 = 58.48° 

n = 1.589 

Z\^ = 180° 

At = 90° 

X = 96.62° 

a = 0.2002 

= 30.48° 

= 59.85° 


3 The Eigenvalue Calibration Method 


The ECM sets out to eliminate/minimise the main sources of error generally 
associated with calibration. Three characteristic matrices are described, W for 
the entrance arm (PSG), A for the polarimeter, and M, the Mueller matrix of 
the sample. The GASP PSG has been modified to that described by|Compain| 


et al. (1999) and has been discussed by Kyne (2014). A measurement of 


I = AMW 


( 2 ) 


is the starting point of the ECM, making use of linear algebra. The ECM 


has the following advantages (Compain et ah, 1999): 


1. No assumption is made to the system to be calibrated - except that it must 
be complete. The precise orientation and position of the various elements 
that comprise the polarimeter do not need to be known for calibration. All 
optical elements are included in the matrix representation. 

2. The characteristics of the reference samples are completely determined dur¬ 
ing the calibration without need for secondary measurements. All aspects 
of the sample are measured during the ECM. W and A, when measured, 
will be a function of the wavelength A, as will the defining characteristics 
of the reference samples. 

3. The accuracy of the calibration procedure is evaluated when the ECM is 
used. 


Therefore, according to the above points, the Stokes vectors for the PSG 
used by GASP does not need to be known prior to calibrating the instrument. 
It is in fact calculated during the calibration process using Equation The 
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following description is based on the work by (Compain et al. 19991 and refer¬ 


ences methods and naming conventions used by (de Martino, 2004 1 and (Lara 


Saucedo 20051. A more detailed description of the mathematical solution for 


this method is found in the above references, including a description of the 
PSG used for calibration in Kyne (2014). 

A linear mapping is described by 


H 4 : M(M) ^ M 4 (K), (3) 

X ^ MX — X{aw)~^{amw). (4) 

M is the Mueller matrix of the reference sample and (amw) corresponds to 
experimental measurements. El has the property of having W within its null 
space, i.e.. 


EI(W) = 0, (5) 

whatever the value of M, because without experimental errors {aw)~^ (amw) 
is equal to W“^MW. A well chosen set of reference samples for {Mi....M„} 
can reduce the number of solutions to Equation to one. And then A can be 
deduced from the following: 


A = {aw)W ( 6 ) 

The chosen reference samples for GASP are based on the work carried out 


by de Martino (2004) and Lara Saucedo (2005) and are discussed in detail by 
these references. 


4 GASP ECM: Methods 

A number of calibration data sets were recorded, while GASP was mounted at 
the Cassegrain focus of the 200 inch Hale telescope at Palomar Observatory. 
An optical lens prescription was used to change the nominal f/16 beam to f/9, 
providing a more acceptable FOV to the GASP PSA. The calibration images 
are acquired as follows: 

— Input light is generated using an LED, a filter, and a linear polariser. 

— An ECM wheel contains a set of samples, AIR, a polariser at 0°, a polariser 
at 90° and a quarter waveplate at 30°. 






















— The polarimeter output intensities are recorded for a set of angles ranging 
from -38° to -209.36° in -10.08° increments for each of the above samples, 
and is repeated using a quarter waveplate (QWP) placed in front of the 
linear polariser. Though of course, as has been explained above, it is not 
necessary to know what these angle are, as long as the same process is 
repeated for the QWP. 

— These intensities are used to extract a system matrix, A and Polarisation 
State Generator (PSG), W. 

The system matrix, A, is normalised to the first element of the theoretical 
GASP PSA using Equationj^to determine variations in the 16 coefficients, and 
a comparison between the individual coefficients. As the experimental matrix 
is calculated based on intensity values, it is difficult to interpret the variation 
in the coefficients without normalisation. 


A 


norm,i 


AiPSA(l,l) 

Ai(l,l) 


(7) 


4.1 GASP EGM: APDs Results 


The EGM is demonstrated in the laboratory using Avalanche Photodiodes, 
where LED light is fed to each detectors, through the polarimeter, by optical 
fibres. A description of the optical setup for these detectors is found in |Kyne| 


(2014). A set of calibration data was recorded, using APDs, in the laboratory 


dated July 2011 and an experimental system matrix, A, shown by, 


A = 


( 0.1766 -0.0901 0.1571 0.0042 \ 

0.4983 -0.2010 -0.4419 0.0113 
0.4992 0.2503 0.1167 -0.4058 
\0.3813 0.1826 -0.0535 0.3229 J 

which is compared to the theoretical PSA, 


( 8 ) 


PSA = 


(9) 


( 0.1766 -0.0858 0.1544 0 \ 

0.1766 -0.0858 -0.1544 0 

0.2176 0.1078 0.0008 -0.1890 
y 0.2176 0.1078 -0.0008 0.1890 ) 

In the absence of errors, and optical configuration differences, these matri¬ 
ces should match. The PSA matrix has been theoretically corrected to produce 
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the system matrix A, which includes polarimetric errors and offsets from the 
design. The matrix is divided into 4 rows. The first 2 rows represent the RPl 
and RP2 channels, and the second 2 rows the TPl and TP2 channels. The 
order of these channels depend on the optical alignment in the laboratory. The 
first column of the system matrix (and the PSA) represents the gain of the 4 
channels of the DOAP, which includes the gain of each detector. The detector 
gain of each APD accounts for the main difference in the overall PSA gain 
coefficients. 

The system matrix A above, was used on a set of data set to reduce the 
Stokes parameters for a polarised light source using an LED, and a glass po- 
lariser (Sodium-Silicate glass). The data was recorded for approximately 48 
hours. The temporal stability of the GASP detectors was not under investiga¬ 
tion. A number of APDs were deemed unsuitable in terms of dark counts and 
electronic noise after these laboratory experiments. 
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Fig. 2 Stokes parameters I, Q, U, and V for an experiment passing LED light through 
a linear polariser to measure the limits of the GASP system using APDs. Each point is 
an average measurement of 1750 seconds of data. The Q, U, and V parameters have been 
normalised to Stokes I to remove any fluctuations in luminosity. A systematic effect is 
observed on Stokes Q and U, the parameters used to the measurement of the degree of 
linear polarisation. A different trend is found for V, which suggests an instrumental effect 
for this laboratory experiment. 


The parameters Q, U, and V have been normalised to Stokes I. A trend in 
I indicates a drop in light intensity, possibly due to an intensity variation in 
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the LED used over the 48 hour period. A change in the laboratory (or LED) 
temperature could also cause this effect, or vibrations due to any movement in 
the laboratory. After normalisation, Q, U, and V show evidence of systematic 
(including random) error over this period. Q, U, and V give mean values of 
0.212 ± 0.003, -0.950 ± 0.002, and 0.015 ± 0.004, respectively. The standard 
deviation value is quoted after normalisation, however, there is a visible trend 
in each of these parameters that suggests a systematic change in the data. This 
is a possibility as during the experiment a linear polariser was used to generate 
a polarised signal. The systematic effect shows a drop in value at time point 
38 (18.5 hours into the experiment) for all Q, U, and V. If the polariser moved, 
and the movement was not symmetric, then this could induce (or reduce) the 
measured linear and circular polarisation. It would also account for the change 
in polarisation angle seen in Figure This systematic error means that the 
errors quoted for this data set are higher than the random noise, which is 
visible in this plot. 

Over the time series the mean degree of linear polarisation (DOLP) was 
measured to be 97.66 ± 0.21%, and the degree of circular polarisation (DOCP) 
0.72 ± 0.36%. This is a promising result that indicates GASP has the ability 
to measure small variations in the linear polarisation error as it fluctuates by 
about 0.2%. The errors quoted are from a combined systematic and random 
effect, however, a polarimetric signal detecting small fluctuations has been 
measured. 

The polarisation angle (PA) was measured prior to the laboratory exper¬ 
iment; this had a value of ^ -38.5°. Figure plots a mean angle of -38.99 ± 
0.19° measured by the ECM with an error level of 0.2°. 

There is variation over the time series ranging over 0.5°. The laboratory 
experiment contains about 48 hours of data in an uncontrolled environment 
subject to temporal fluctuations and potential mechanical disturbances. The 
temperature was found to be variable over the course of the experiment, but 
this data was not recorded. It is also noted that structure holding the linear 
polariser was subject to some movement within the mounting device, but it 
is not possible to prove that this is what caused the change in PA. If the 
variation in the PA is examined on a shorter time scale - in an area where the 
systematic effect is reduced, or absent - the variation is less than 0.1°. This is 
a good indication that GASP is working correctly, and it also illustrates the 
capability of GASP to detector fluctuations in the PA to a level of ± 0.1°. 
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Fig. 3 Polarisation angle measurement for an experiment passing LED light through a 
linear polariser to measure the limits of the GASP system using APDs. The drop, and 
variation, in Stokes Q and U is observed in this plot. There is a variation of 0.5° from 
beginning to end, which suggests a number of errors including thermal or pressure changes, 
or that the polariser moved during the experiment. This is not a random effect and over a 
prolonged period (2 days) it is possible that there was movement in the laboratory, or in 
the vicinity of the instrument to account for this variation. 


4.2 GASP ECM: EMCCDs Results 


There were obstacles encountered using the APDs from Section 4.1 which 
included failing modules and high noise levels and dark counts. 2 Andor Ultra 
897 EMCCD detectors were loaded to the project by a collaborating group at 
the California Institute of Technology, therefore, GASP was modified for use 
as an imaging polarimeter. 

A number of observation targets were used to produce a geometrical trans¬ 
formation of the RPl channel with each of the RP2, TPl, and TP2 channels 


(Kyne 2014). Tyo et al. (2006) and Smith et al. (1999) discuss the critical issue 


of image registration for sequential or simultaneous image collection. Regis¬ 
tration to 1/20 of a pixel can achieve less than 0.01 error in DOLP and DOCP 
results. Misregistration can result, for example, from separate focal planes that 
are not looking at the same region of space; it can also result in beam wan¬ 
der from a rotating element. In practice, achieving even a half-pixel alignment 
mechanically is difficult - of course this will be dependent on the pixel size of 
the detector - and software post-processing alignment is frequently necessary. 
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GASP overcomes the problem of motion of the target held across the FOV 
using the DOAP technique. Multiple images per channel are acquired simul¬ 
taneously, removing the need for continuous guiding. However, spatial reg¬ 
istration of multiple images is complicated by the need to correct for both 
mechanical misalignment as well as optical ‘misalignment’ arising from aber¬ 
rations due to separate optical paths. The TP has a longer optical path length 
compared to RP, and it was found in Collins et al. (2009) that a design feature 
of the transmitted path resulted in an astigmatism after transmission and be¬ 
fore reimaging. It was not possible to correct this aberration for the purpose 
of these observations. 

Registration is carried out using matched points, channel-to-channel, from 
a set of targets on Night 4 of the Palomar observing run. Kyne (2014) gives a 
detailed explanation for the point selection process using observational data. 
The registration routine involved the following steps: 


1. A reference channel is selected - RPl is chosen for this data analysis. 

2. A set of points are matched in a region of interest (ROI) to be calibrated 
(analysed) between RPl and the remaining channels. 

3. geomap, an IRAF command, is used to generate a geometric transformation 
between RPl and the remaining channels based on this matching. 

4. gregister (from IRAF) is then used to register the images. 

The central region of the images in Figure [^appear to be well registered, 
however, looking from the centre out to the edges of the field the matching 
deteriorates. This is expected given the limited number of matched points, 
and overall optical image quality. A pixel-by-pixel polarimetric analysis was 
carried out on an area of the FOV that showed minimum distortion. 

Figurej^is an image of the normalised calibration system matrix generated 
from this registration solution. Each pixel has its own coefficient which forms 
a separate system matrix, for all pixels in this area of the FOV. In this way the 
analysis can be compared to that used for APDs, as each pixel is be thought 
of as a single APD. The main difference in this calibration analysis is that the 
result is field dependent, as when the pupil is reimaged it is spread across a 
number of pixels, instead of a single APD. 

It has been discussed that the hrst column of the matrix gives the gain for 
each channel. The first row is a measure for RPl, row 2-4, are RP2, TPl, and 
TP2, respectively. There is a variation (banding) in the gain values, showing a 
systematic trend across each coefficient in column 1, with the exception of gain 
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Fig. 4 All registered GASP channels using science data points and image artefacts. A 
pixel-by-pixel calibration analysis will be analysed on the area indicated by the orange box. 
The shape of the RP2, TPl, and TP2 images are not perfectly matched to RPl due to 
insufficient matched points. IT is also noted that the optical aberration on the TP is too 
severe to correct in software alone. Top LHS: RPl, top RHS: RP2, bottom LHS: TPl and 
bottom RHS: TP2. 


coefficient RPl. On examination of the trends found in the image residuals 
after registration (Kyne, 20141, the patterns observed are nearly identical. 
It was not possible to completely eliminate all optical distortion using post¬ 
processing techniques, resulting in these systematic effects. Normalisation gives 
a good indication of how well registered the RP2, TPl, and TP2 channels are 
compared to RPl. The system matrix is used to perform the Stokes reduction 
for science data in Section [5l 

The gain values for each pixel for RPl and RP2 compare well, as do those 
for TPl and TP2 getting worse towards the edge of the field. These values 
are compared to the theoretical PSA, though, it is not absolutely correct to 
compare this to the theoretical PSA given the differences that will exist in 
optical alignment. In the absence of variation in luminosity in the field, of 
each channel, a more random pattern would be observed in each coefficient. 
Where the values are expected to be low (~ 0), a pattern will be less evident 
where the random errors are approximately equal to the systematic, taking 
into account alignment variation in the theoretical PSA. 
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Fig. 5 A pixel-by-pixel experimental calibration matrix for the image section (~ 10 ") 
indicated by the orange box in Figure]^ The calibration images were registered to channel 
RPl using science data. 


4.3 GASP ECM: Verification 

It is important that the results of the ECM are verified experimentally to 
ensure that the system is calibrated. This is done by testing the system matrix 
with another set of calibration data, which is passed through the same samples 
used for the ECM. A number of parameters are checked in this way: the 
polarimetric offsets/error as a function of pixel, or subsection of the FOV, and 
the field-dependence. Each ECM sample has been aligned to a particular value, 
which the ECM measures as an output of the calibration procedure. There is 
optical field distortion present in the system, which could not be completely 
removed in software; we can map the calibration, pixel-to-pixel variation over 
the FOV by this verification process. 

A systematic trend can be found across the images of the coefficients in 
Figure which correlates with the patterns visible in the raw data after 
registration. The Mueller matrix for the AIR ECM sample is found in Figure 
1^ The a values are as low as 0.5%, compared with the expected MM. There 
is evidence of a systematic effect on each of the TP images, particularly TP2. 
This is consistent with the fact that there is more observed distortion on this 
path. The theoretical and mean value for the second to last coefficient is quite 
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high; this is indicative of the difficulty in image registration. The same is found 
for the last TPl coefficient. 
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Fig. 6 The normalised Mueller matrix for the AIR sample for the image section indicated 
in Figure]^ The calibration images were registered to channel RPl using science data. 


Table 2 Theoretical and mean values for the normalised Mueller matrix of the AIR ECM 
sample; registration performed using science data points. 


mm air 


X 


/1.0000 0.0000 0.0000 0.0000 \ 

0.0000 1.0000 0.0000 0.0000 

0.0000 0.0000 1.0000 0.0000 

\ 0.0000 0.0000 0.0000 1.0000 / 


/ 1.0000 0.0004 0.0014 -0.0020 \ 

0.0043 1.0021 -0.0018 -0.0045 

-0.0032 -0.0008 0.9854 0.0467 

\-0.0094 0.0046 0.0873 0.9902 / 


Figure and are the Mueller matrices for the linear polariser samples 
PO and P90, respectively. These samples have similar MM, with the exception 
of sign. The variation is low on each channel, however, the values for the 
3'’‘^ coefficients on rows 1 and 2 are higher than expected. This could be for a 
number of reasons, image registration, imperfections in the polarisers used, the 
polarisers may not exactly align to the stated value. The pixel values across 
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the FOV will change as a function of angle of incidence between the sample 
and the detector - which is related to the quality of the reimaging optics. 

There is a trend on each of the TP rows (3 and 4). The top left and bottom 
right corner show larger variation compared to that in the centre of the FOV. 
It is most likely due to an absence of registration points in those areas. The 
PO and P90 experimentally determined Mueller matrices show low variation 
with their corresponding theoretical matrices. 
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Fig. 7 The normalised Mueller matrix for the PO sample for the image section indicated in 
Figure The calibration images were registered to channel RPl using science data. 


Table 3 Theoretical and mean values for the normalised Mueller matrix of the PO ECM 
sample; registration performed using science data points. 


IVIlVIpQ 


X 


/ 1.0000 1.0000 0.0000 0.0000 \ 

1.0000 1.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 

\ 0.0000 0.0000 0.0000 0.0000 / 


/ 1.0000 1.0075 -0.0166 0.0015 \ 

0.9848 1.0169 -0.0175 0.0017 

-0.0004 -0.0014 -0.0011 0.0015 

\ 0.0044 0.0049 -0.0000 0.0003/ 


The results found for the horizontal polariser in Figure are comparable 
to those in Figure and the mean coefficient values can be found in Table 
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and 1^ respectively. The variation is similar and the proposed systematic 
effect of registration is also present in each of these images. These Mueller 
matrices show a variability compared to the field-dependence errors found in 
the calibration matrix (Figure]^. 
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Fig. 8 The normalised Mueller matrix for the P90 sample for the image section indicated 
in Figure]^ The calibration images were registered to channel RPl using science data. 


Table 4 Theoretical and mean values for the normalised Mueller matrix of the P90 ECM 
sample; registration performed using science data points. 


MMpgo 

X 


/ 1.0000 - 1.0000 0.0000 0.0000 \ 

- 1.0000 1.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 j 


/ 1.0000 - 0.9930 0.0238 - 0.0002 \ 
- 1.0091 0.9969 - 0.0242 - 0.0005 
0.0190 - 0.0181 0.0012 - 0.0017 
1 ^- 0.0017 0.0010 - 0.0002 - 0 . 0002 ^ 



Figure is the Mueller matrix for the sample R30. It was concluded, from 
an analysis of averaged/binned data, that this sample has been measured to 
be closer to an angle of ~ 38°. The sample was originally aligned to approx¬ 
imately ~ 30°, a difficult mechanical task with the current setup, however, 
when calibration was carried out for a point source (a pinhole at the instru- 
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ment focus, instead of an extended field) the ECM measured this sample to 
be ~ 38°, which is another very useful attribute of the method. Just like the 
figures above, there is a low systematic trend present which can be compared 
to that in the system matrix. 

Looking at Figure the TP coefficients, rows 3 and 4, register the highest 
variation for all coefficients, about 1%, which compares to the FOV variation 
for that of the system matrix (Figure]^. As mentioned above, this is most 
likely a result of incorrect pixel-matching. 
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Fig. 9 The normalised Mueller matrix for the R30 sample for the image section indicated 
in Figure]^ The calibration images were registered to channel RPl using science data. 


Table 5 Theoretical and mean values for the normalised Mueller matrix of the R30 ECM 
sample; registration performed using science data points. 


]VI]VIp 38 X 


/ 1.0000 0.0000 0.0000 0.0000 \ 
0.0000 0.0585 0.2347 -0.9703 
0.0000 0.2347 0.9415 0.2419 

0.0000 0.9703 -0.2419 0.0000 ^ 


/ 1.0000 -0.0010 0.0023 -0.0014 \ 
0.0083 0.0212 0.2259 -0.9689 

-0.0040 0.2593 0.9217 0.2823 

-0.0020 0.9720 -0.1346 -0.0545^ 
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The variation from the mean value for each coefficient of each Mueller 
matrix is low, and within the expected values for each of the RPl, RP2, TPl, 
and TP2 channels, which was discussed in Section |4.2[ A pattern remains in 
the FOV of a number of coefficients in each of the samples for AIR, PO, P90, 
and R30. The same pattern is found in Figure for the respective channels. 
The reason for this pattern can also be referenced in Figure where there are 
insufficient points, and the areas absent of points are showing this trend. 


5 Zenith Flat-Field Observations 


A 35 minute set of data was recorded at twilight on the last night of observa¬ 
tions during the Palomar November observing run. Each frame was integrated 
for 1 second, continuously, as the Sun was rising. The data was recorded in 
this way to measure how the polarisation of the sky changes with respect to 
the elevation of the Sun over time, i.e. the polarisation of the sky is expected 
to increase as the elevation increases. This type of experiment gives valuable 
information about the accuracy of the GASP image registration/calibration 
by measuring the limit of the polarimetric variation over the FOV. It is also a 
very interesting experiment to determine a measurement of the degree of linear 
polarisation, and change in polarisation angle, of the twilit sky over time. 


Dahlberg et al. (2009) and Harrington et al. (2011) describe how there are 


many atmospheric and geometric considerations, which determine the skylight 
polarisation at a particular observation site. The magnitude of the polarisa¬ 
tion angle can depend on the solar elevation, atmospheric aerosol content, 
aerosol vertical distribution, aerosol scattering phase function, wavelength of 
the observation and secondary sources of illumination. Further descriptions of 


this can be found in these references Horvath et al. (2002); Suhai and Horvath 


(2004); Cronin et al. (2006). Anisotropic scattered sunlight that can arise from 


reflections off land or water can be highly polarised and variable |PeltoniemI| 


et al. (2009); Litvinov et al. (2010); Kisselev and Bulgarelli (2004). Aerosol 


particle optical properties and vertical distributions in the atmosphere can 


vary and can cause variability on the incoming polarisation signal (Shukurov 


[and Shukurov[ [2^006[ [Vermeulen et al.[ |2000] [Ugolnikov et al.[|2004[[Ugolnikov| 


and Maslov, 2010). The polarisation can change across atmospheric absorp¬ 


tion bands or can be influenced by other scattering mechanisms Boesche et al. 


(2006); Zeng et al. (2008); Aben et al. (1999, 2001). The polarimetric mea¬ 


surements made by GASP (Harrington et al. 2011), will contain deviations 
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from the Rayleigh Sky model, which will grow with aerosol, cloud, ground or 
sea-surface scattering, and affect the telescope line-of-sight. However, clear, 
cloudless, low-aerosol (ideal) conditions should yield high linear polarisation 


amplitudes and small deviations in polarisation, from the model (Pust and 


Shaw 

2005 

2006 

2007 

2008 2009 

Shaw et al. 

2010 


This analysis will assess the performance of the GASP calibration, as well 
as carrying out scientific measurements. The Zenith data was recorded in 3 
stages: The first was summed and averaged for 48 seconds, the first image in 
the figures to follow. Images 2-11 are averaged over 200 frames (3.3 minutes). 
The final image, the brightest sky image, is 134 seconds of data. 


5.1 Zenith Flat-Field Observations: Polarimetric Results 


An approximation for the expected variation in polarisation angle related to 
time-per-frame is give in Table The variation in DOLP is compared with 
that of the Rayleigh Sky Model, which gives the value for the degree of linear 
polarisation calculated as a function of the absolute angle for a reflected beam 
of light, 


D0LP(7) = 


1 — cos^ 7 


( 10 ) 


1 -I- cos^ 7 

7 is the scattered angle formed between the (telescope) pointing and the 
Sun, and 9 = 90° — 7 , the angle made by the Sun, the pointing and the 
Zenith. According to Equation [TOl as the scattered angle increases, the DOLP 
also increases. A detailed experiment of how this model describes the E-vector 
was carried out by Suhai and Horvath (2004). 


Table 6 Expected polarisation angle variation over the same 35 minute period of observa¬ 
tion by GASP for Zenith flat-field data. The first time is not usable as we are measuring 
from this value; The APA values are measured from here. 


Flat-field 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Time (mins) 

0.47 

2.45 

5.75 

9.05 

12.35 

15.65 

18.95 

22.25 

25.55 

28.85 

32.15 

35.55 

ZiPA (°) 

0.42 

0.69 

0.69 

0.69 

0.69 

0.69 

0.69 

0.69 

0.69 

0.69 

0.71 


The values for time in Table were calculated from the time-stamps in 
the header (data information) files. The Sun moves 1° in 4 minutes and, given 
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that the elevation at Palomar is +32.3°, the expected change in PA is 0.84° in 
4 minutes. The times given in Table are the midpoints of the recorded data 
- taking into account delays in acquisition due to stoppages. 


5.2 Discussion 


The results for the DOLP are found in Figure Some of the patterns visible 
in the system matrix, in Figure appear in the same location for the gain 
coefficients on RP2, TPl, and TP2. Overall, a flat region of linear polarisation 
of this section of sky is measured with a minor gradient from the bottom-left- 
corner to the top-right-corner. The largest value for a per frame is about 1% 
for frames with higher SNR. 

The GASP system has been compared to work carried out by [Smith et al.| 
(1999) where, if GASP measures a pixel error of ~ 0.3 pixels (propagation of 


the RMS errors in all 3 channels for 4 x 4 binning), this could result in ± 
6% error in the DOLP and DOGP. An error for the PA has not been stated 
by the authors. There is variation in the pixel-to-pixel values measured by the 
GASP system that are comparable to this result in the extended FOV, outside 
of that in Figure]^ 


22 






1. mean = 35.63% 


2. mean = 35.17% 


3. mean = 31.71% 



2468 10 2468 10 2468 10 


Fig. 10 The degree of linear polarisation as a function of time, and field, are plotted. 
These images are produced using a registration method by science data points. The DOLP 
is measured as the Sun rises; the first was summed and averaged for 48 seconds, images 2 - 
11 are averaged over 200 frames (3.3 minutes), and the final image, the brightest sky image, 
is 134 seconds of data. There is a clear change in pattern as the sky brightens, and overcomes 
the noise floor. The mean and a values are given in the title of each flat. The image has been 
flipped and rotated compared to the images in Figure to match the telescope correction. 
This is the reason for the blue pixels at each corner of the image. The PA vectors are also 
plotted. The cr value refers to the random and systematic error over the FOV of interest in 
Figure]^ for the 3x3 pixel area indicated by the rectangle on the final image. 
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1. mean= 88.92° 



2. mean= 84.29° 




7. rr^afi = 4l.3'P° 8. n^afi = 4^.39° 9. rt^afi = 31.39° 
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Fig. 11 The polarisation angle as a function of time, and field, are plotted. These images 
are produced using a registration method by science data points. The PA is measured as 
the Sun rises; the first was summed and averaged for 48 seconds, images 2-11 are averaged 
over 200 frames (3.3 minutes), and the final image, the brightest sky image, is 134 seconds 
of data. There is a clear change in pattern as the sky brightens, and overcomes the noise 
floor. The mean and a values are given in the title of each flat. The image has been flipped 
and rotated compared to the images in Figure]^ to match the telescope correction. This is 
the reason for the blue pixels at each corner of the image. The cr value refers to the random 
and systematic error over the FOV of interest in Figure]^ for the 3x3 pixel area indicated 
by the rectangle on the final image. 
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Fig. 12 The degree of circular polarisation as a function of time, and field, are plotted. 
These images are produced using a registration method by science data points. The DOCP 
is measured as the Sun rises; the first was summed and averaged for 48 seconds, images 2 - 
11 are averaged over 200 frames (3.3 minutes), and the final image, the brightest sky image, 
is 134 seconds of data. There is a clear change in pattern as the sky brightens, and overcomes 
the noise fioor. The mean and a values are given in the title of each flat. The image has been 
flipped and rotated compared to the images in Figure to match the telescope correction. 
The a value refers to the random and systematic error over the FOV of interest in Figure]^ 
for the 3x3 pixel area indicated by the rectangle on the final image. 
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A similar trend in the pattern for the polarisation angle (PA) is found in 
Figure |11[ The mean and a values are taken from the 3x3 pixel area high¬ 
lighted in the last frame in Figure [lOl A low error value has been found in a 
reasonably flat section of the field, comparable to the same area in the polari- 
metric measurements for the DOLP and DOCP. The same areas of the FOV as 
those for DOLP show a gradual gradient in the pattern, particularly from the 
bottom left to top right corner. This area is also present in the registered TP2 
images and appears in the calibration matrix. It is noted that the polarimetric 
results have been flipped and rotated to account for the telescope orientation 
on the sky. This has not been performed for the calibration measurements. 
The variation in pattern is a close match to the image of the TPl gain coeffi¬ 
cient in Figure A comparison is also made between the gain coefficient for 
the polarisation angle and RP2. This makes sense as the gain coefficient for 
RP2 and TPl are used to calculate the values for Stokes Q and U respectively 
when calculating the PA. 

The first 6 frames show a combination of random and systematic error in 
the results for the DOLP. This has not changed when measuring the PA and 
the conclusions are the same in terms of the calibrated gain of the system. 
Table 1^ shows that over a 35 minute period the expected change in angle of 
the Sun with respect to the Earth is ~ 6.3°. 

The pattern in the images for the DOCP in Figure [T^ also show systematic 
variation. The highlighted area in the final frame measures the DOCP at 0.7 ± 
0.3% and appears relatively flat, which indicates that GASP is measuring a low 
level of circular polarisation, almost zero, and a reasonable a variation. This 
area is used as outside the total blue section there is an indication that the field 
is less flat (where the colour map changes from blue to red, left-to-right). There 
is a large increase in DOCP and the visible shape in the bottom right corner 
is altered compared to the left-hand-side. This corner is showing quite a large 
variation in optical distortion not corrected by geometrical registration. It is 
also noted that this pattern exactly matches that of the TP2 gain coefficient in 
Figure which measures the Stokes V parameter. This pattern is also found 
in the raw image and flat-held residuals for TP2 images. 

The results in Tablej^are measurements of the mean in the highlighted 3 x 
3 pixel area in the final image of each polarimetric measurement. A variation 
in the polarisation angle of ^ 52° over the course of the observation. Realis¬ 
tically, this value should be measured for when there are less influences from 
random variations as a result of detector noise. These initial measurements of 
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Table 7 Mean polarimetric results for Zenith fiat-field data using a pixel-by-pixel analysis 
for the 3x3 pixel area in the final image of each polarimetric measurement. The APA 
values are also given. 


Flat-field 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Time (mins) 

0.47 

2.45 

5.75 

9.05 

12.35 

15.65 

18.95 

22.25 

25.55 

28.85 

32.15 

35.55 

DOLP (%) 

35.63 

35.17 

31.71 

35.67 

44.36 

57.40 

69.63 

79.23 

84.10 

86.01 

87.93 

90.02 

DOCP (%) 

1.62 

1.94 

1.19 

0.41 

-0.49 

-1.32 

-1.84 

-1.97 

-1.91 

-1.69 

-1.44 

-1.23 

PA(“) 

88.92 

84.29 

73.23 

62.55 

52.55 

45.67 

42.57 

40.61 

39.25 

38.09 

37.07 

36.24 

APA (“) 

- 

4.63 

11.06 

10.68 

10.00 

6.88 

3.10 

1.96 

1.36 

1.16 

1.02 

0.83 


DOLP are mostly at the noise level, which will cause the polarisation to be 
overestimated. There appears to be a significant increase in the light level at 
frames 5-6, which gives a change in the PA of ^ 7°. The lowest frame-to- 
frame variation in PA observed is 0.83 ± 0. lOf] The APA value is decreasing 
over time, approaching 0.83°, compared with an expected variation of 0.71°. 
It is possible that if more data was recorded the APA would terminate at this 
value. Evidence for a dependence on a high level of counts per pixel is clear, 
however, more data is required over a longer observing period. 

The value for the DOLP increases in value between frames 4-9, where 
there is a significant change in light level. However, the A value, frame-to- 
frame, for the DOLP reaches a terminal value between frames 9 and 10. From 
frame 9 onwards, the ADOLP is 2%, though with limited data it is difficult to 
comment on this trend. The change in DOLP over time, and as a function of 
PA shows good agreement with the Rayleigh Sky Model, however, the experi¬ 
mental measurements indicate that something in addition to Rayleigh scatter 
is detected by GASP. The measured DOCP value changes by 1 - 2% around 
the expected value of 0%, eventually converging at 0%. This change is slow, 
by 0.2% in the hnal frames reaching a value of —1.23 ± 0.3%. Deviations from 
the Rayleigh Sky model have already been discussed, and given the nature 
of these observations, it is unlikely that GASP has measured pure Rayleigh 
scatter given the poor observing conditions which include poor seeing, cloud 
cover and unknown scatter from other local sources or contaminants. However, 
in all cases each result indicates that GASP is measuring a field dependent, 
changing polarimetric signal. 


^ Errors quoted are based on the Poissonian statistics for the 3x3 pixel area 
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6 Conclusion 


The use of the ECM for an APD detector configuration indicate that GASP 
works well as a division of amplitude polarimeter. This was concluded as a 
comparable system matrix to the theoretical PSA was obtained and tested on 
linearly polarised light. The ECM measured absolute values for the DOLP, 
DOCP, and PA, with variation in these values as low as 0.2% and 0.1° (in 
the absence of systematic error), respectively. This experiment also showed a 
systematic error in the Stokes parameters, which indicates that GASP was 
measuring a possible mechanical movement, or change in temperature, in the 
system. 

Calibrating GASP using EMCCDs was a more difficult process, given the 
optical system that was needed to adapt the instrument for observing on the 
200 inch at Palomar. The main difficultly was due to image registration when 
post-processing the data, however, a simple laboratory experiment can verify 
that these errors are low (< %1), but optical distortion (after image registra¬ 
tion) remains at the edge of the field. These verification results show that it 
was possible to calibrate GASP, with error; the calibration was further tested 
on a set of data recorded while the telescope was pointed at the Zenith. It 
was found that there was little variability in the spatial polarimetric signal 
across the FOV and that, in spite of registration errors, it was possible to 
reduce a value for the polarisation angle of the Sun based on its geometric po¬ 
sition relative to the Earth. These results are comparable to those predicted by 
the Rayleigh Sky Model, however, further observations, for extended periods, 
and higher SNR, are required to investigate/verify the results found from this 
experiment. 


References 

Aben, L, Helderman, F., Stam, D. M., and Stammes, P. (1999). Spectral 
fine-structure in the polarisation of skylight. Geophysical Research Letters, 
26:591-594. 

Aben, L, Stam, D. M., and Helderman, F. (2001). The ring effect in skylight 
polarisation. Geophysical Research Letters, 28:519-522. 

Azzam, R. and Bashara, N. (1988). Ellipsometry and Polarized Light. North 
Holland. 


28 



Azzam, R. M. A. and Lopez, A. G. (1989). Accurate calibration of the four- 
detector photopolarimeter with imperfect polarizing optical elements. Jour¬ 
nal of the Optical Society of America A, 6(10):1513-1521. 

Boesche, E., Stammes, P., Ruhtz, T., Preusker, R., and Fischer, J. (2006). 
Effect of aerosol microphysical properties on polarization of skylight: sensi¬ 
tivity study and measurements. Applied Optics, 45:8790-8805. 

Collins, P., Kyne, G., Lara, D., Redfern, M., Shearer, A., and Sheehan, B. 
(2013). The Galway astronomical Stokes polarimeter: an all-Stokes optical 
polarimeter with ultra-high time resolution. Experimental Astronomy. 

Collins, P., Sheehan, B., Redfern, M., and Shearer, A. (2009). GASP - Galway 
Astronomical Stokes Polarimeter. ArXiv e-prints. 

Compain, E. and Drevillon, B. (1998). Broadband Division-of-Amplitude Po¬ 
larimeter Based on Uncoated Prisms. Applied Optics, 37:5938-5944. 

Compain, E., Poirier, S., and Drevillon, B. (1999). General and Self-Consistent 
Method for the Calibration of Polarization Modulators, Polarimeters, and 
Mueller-Matrix Ellipsometers. Applied Optics, 38:3490-3502. 

Cronin, T. W., Warrant, E. J., and Greiner, B. (2006). Celestial polarization 
patterns during twilight. Applied Optics, 45:5582-5589. 

Dahlberg, A. R., Pust, N. J., and Shaw, J. A. (2009). All-sky imaging po¬ 
larimeter measurements of visible and NIR skylight at Mauna Loa, Hawaii. 
In Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, volume 7461 of Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, page 7. 

de Martino, A. (2004). General methods for optimized design and calibration 
of Mueller polarimeters. Thin Solid Films, 455:112-119. 

Drevillon, B. (1993). Progress in crystal growth and characterization of ma¬ 
terials. 

Harrington, D. M., Kuhn, J. R., and Hall, S. (2011). Deriving Telescope 
Mueller Matrices Using Daytime Sky Polarization Observations. Publica¬ 
tions of the Astronomical Society of the Pacific, 123:799-811. 

Hauge, P. S. (1978). Mueller matrix ellipsometry with imperfect compensators. 
Journal of the Optical Society of America A, 68(11):1519-1528. 

Horvath, G., Barta, A., Gal, J., Suhai, B., and Haiman, O. (2002). Ground- 
based full-sky imaging polarimetry of rapidly changing skies and its use for 
polarimetric cloud detection. Applied Optics, 41:543-559. 

Kisselev, V. and Bulgarelli, B. (2004). Reflection of light from a rough water 
surface in numerical methods for solving the radiative transfer equation. 


29 



Journal of Quantitative Spectroscopy and Radiative Transfer, 85:419-435. 

Kyne, G. (2014). The optical development and ealibration of the Galway As¬ 
tronomical Stokes Polarimeter (GASP) as a multi-detector system for the 
polarimetrie observations of variable optical sources. PhD thesis, The Na¬ 
tional University of Ireland, Galway, http://hdl.handle.net/10379/4572. 

Lara Saucedo, D. (2005). Three-dimensional Gomplete Polarisation Sensitive 
Imaging using a Gonfocal Mueller Matrix Polarimeter. PhD thesis. Imperial 
College London. 

Litvinov, P., Hasekamp, O., Cairns, B., and Mishchenko, M. (2010). Reflec¬ 
tion models for soil and vegetation surfaces from multiple-viewing angle 
photopolarimetric measurements. Journal of Quantitative Speetroscopy and 
Radiative Transfer, 111:529-539. 

Peltoniemi, J., Hakala, T., Suomalainen, J., and Puttonen, E. (2009). Polarised 
bidirectional reflectance factor measurements from soil, stones, and snow. 
Journal of Quantitative Spectroseopy and Radiative Transfer, 110:1940- 
1953. 

Pust, N. J. and Shaw, J. A. (2005). Dual-field imaging polarimeter for studying 
the effect of clouds on sky and target polarization. In Shaw, J. A. and Tyo, 
J. S., editors. Polarization Seience and Remote Sensing II, volume 5888 
of Society of Photo-Optical Instrumentation Engineers (SPIE) Conferenee 
Series, pages 295-303. 

Pust, N. J. and Shaw, J. A. (2006). Dual-field imaging polarimeter using liquid 
crystal variable retarders. Applied Optics, 45:5470-5478. 

Pust, N. J. and Shaw, J. A. (2007). All-sky polarization imaging. In Soci¬ 
ety of Photo-Optieal Instrumentation Engineers (SPIE) Oonference Series, 
volume 6682 of Society of Photo-Optieal Instrumentation Engineers (SPIE) 
Oonference Series, page 4. 

Pust, N. J. and Shaw, J. A. (2008). Digital all-sky polarization imaging of 
partly cloudy skies. Applied Optics, 47:190. 

Pust, N. J. and Shaw, J. A. (2009). How good is a single-scattering model of 
visible-NIR atmospheric skylight polarization? In Society of Photo-Optical 
Instrumentation Engineers (SPIE) Oonference Series, volume 7461 of Soci¬ 
ety of Photo-Optieal Instrumentation Engineers (SPIE) Oonference Series, 
page 0. 

Shaw, J. A., Pust, N. J., Staal, B., Johnson, J., and Dahlberg, A. R. (2010). 
Continuous outdoor operation of an all-sky polarization imager. In Soci¬ 
ety of Photo-Optieal Instrumentation Engineers (SPIE) Oonference Series, 


30 



volume 7672 of Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, page 0. 

Shukurov, A. K. and Shukurov, K. A. (2006). Field studies of the correlation 
between the atmospheric aerosol content and the light polarization at the 
zenith of the daytime sky. Izvestiya Atmospheric and Oceanic Physics, 
42:68-73. 

Smith, M. H., Woodruff, J. B., and Howe, J. D. (1999). Beam wander consid¬ 
erations in imaging polarimetry. In Goldstein, D. H. and Chenault, D. B., 
editors. Society of Photo-Optical Instrumentation Engineers (SPIE) Con¬ 
ference Series, volume 3754 of Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, pages 50-54. 

Suhai, B. and Horvath, G. (2004). How well does the Rayleigh model de¬ 
scribe the E-vector distribution of skylight in clear and cloudy conditions? 
A full-sky polarimetric study. Journal of the Optical Society of America A, 
21(9):1669-1676. 

Thompson, R. C., Bottiger, J. R., and Fry, E. S. (1980). Measurement of po¬ 
larized light interactions via the mueller matrix. Applied Optics, 19(8):1323- 
1332. 

Tyo, J. S., Goldstein, D. L., Ghenault, D. B., and Shaw, J. A. (2006). Review 
of passive imaging polarimetry for remote sensing applications. Appl. Opt., 
45(22):5453-5469. 

Ugolnikov, O. and Maslov, I. (2010). Quadrantids 2008 and 2009: Detection 
of Dust in the Atmosphere by Polarization Twilight Sky Measurements. In 
Proceedings of the International Meteor Conference, 27th IMC, Sachticka, 
Slovakia, 2008, pages 98-101. 

Ugolnikov, O. S., Postylyakov, O. V., and Maslov, I. A. (2004). Effects of mul¬ 
tiple scattering and atmospheric aerosol on the polarization of the twilight 
sky. Journal of Quantitative Spectroscopy and Radiative Transfer, 88:233- 
241. 

Vermeulen, A., Devaux, C., and Herman, M. (2000). Retrieval of the Scatter¬ 
ing and Microphysical Properties of Aerosols from Ground-Based Optical 
Measurements Including Polarization. I. Method. Applied Optics, 39:6207- 
6220. 

Zeng, J., Han, Q., and Wang, J. (2008). High-spectral resolution simulation of 
polarization of skylight: Sensitivity to aerosol vertical profile. Ceophysical 
Research Letters, 35:20801. 


31 



